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The dynamics of monopole giant resonances in nuclei is analyzed in the 
time-dependent relativistic mean-field model. The phase spaces of isoscalar 
and isovector collective oscillations are reconstructed from the time-series of 
dynamical variables that characterize the proton and neutron density distri- 
butions. The analysis of the resulting recurrence plots and correlation dimen- 
sions indicate regular motion for the isoscalar mode, and chaotic dynamics 
for the isovector oscillations. Information-theoretic functionals identify and 
quantify the nonlinear dynamics of giant resonances in quantum systems that 
have spatial as well as temporal structure. 

I. INTRODUCTION 

In a recent article |]J we have started the analysis of collective nonlinear dynamics in 
atomic nuclei in the framework of time-dependent relativistic mean field theory. Atomic 
nuclei provide excellent examples of quantum systems in which the transition from regular 
to chaotic dynamics can be observed in experiments, and studied with a variety of sophisti- 
cated theoretical models. Signatures of chaotic dynamics have been observed in correlations 
of nuclear level distributions, and in the microscopic and collective motion of the nuclear 
many-body system ||. Especially interesting in this respect are giant resonances: highly 
collective nuclear excitations whose properties, excitation energies and widths, nevertheless 
reflect the underlying microscopic dynamics. Theoretical studies predict that regular col- 
lective modes coexists with chaotic single-nucleon motion: the adiabatic mean-field created 
by the nucleons averages out the random components of their motion. In Ref. |]J we have 
studied the dynamics of the most simple giant resonances: isoscalar and isovector collective 



monopole oscillations. In these resonances only spatial degrees of freedom are excited and 
the motion is spherically symmetric, and therefore relatively simple for numerical integra- 
tion. We have analyzed time-dependent and self-consistent calculations that reproduce the 
experimental data on monopole giant resonances in spherical nuclei. In the microscopic 
mean-field model, self-consistent solutions for ground-states provide initial conditions, and 
fully time-dependent calculations are performed for the single-nucleon dynamics. Due to the 
self-consistent time-evolution, the nuclear model system is intrinsically non-linear. In par- 
ticular, we have studied the difference in the dynamics of isoscalar and isovector collective 
modes. Time-series, Fourier power spectra, Poincare sections, autocorrelation functions, and 
Lyapunov exponents have been used to characterize the nonlinear system and to identify 
chaotic oscillations. It has been shown that the oscillations of the collective coordinate can 
be characterized as regular for the isoscalar mode, and that they become chaotic when initial 
conditions correspond to the isovector mode. 

In the present article we continue our investigation of giant monopole resonances. There 
are two main objectives. Firstly, the analysis of phase spaces reconstructed from time-series 
of collective dynamical variables that characterize the isoscalar and isovector oscillations. 
The structure of recurrence plots and the correlation dimensions should provide additional 
information about the chaotic regime of collective motion. The second objective is to study 
possible applications of modern information-theoretic techniques to the nuclear many-body 
dynamics (analysis of information entropies and evaluation of mutual information functions). 
This is especially important since nuclei are quantum objects of finite size, and therefore 
display interesting spatial as well as temporal behavior. We will discuss several information 
functionals that could be used to examine the influence of the finite spatial extension of 
nucleon densities on the nonlinear dynamics of collective excitations. 

The article is organized as follows. In Sec. 2 we describe the time-dependent relativistic 
mean-field model. The time-series analysis and the reconstruction of phase spaces are per- 
formed in Sec. 3. The use of information-theoretic functionals in the description of nonlinear 
dynamics is discussed in Sec. 4. A summary of our results is presented in Sec. 5. 
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II. TIME-DEPENDENT MEAN-FIELD MODEL 



Theoretical models based on quantum hadrodynamics ||,U have been remarkably suc- 
cessful in the description of many physical phenomena in atomic nuclei, nuclear matter, 
neutron stars, heavy-ion collisions, electron scattering on nuclei. In particular, the relativis- 
ts mean-field model has been applied in calculations of properties of ground and excited 
states, both for spherical and deformed nuclei (for a recent review see ||). In Refs. @-|^] 
we have used the time-dependent version of the relativistic mean-field model to describe 
the dynamics of giant resonances. These collective excitations have been observed in nuclei 
over the whole periodic table, and their characteristic properties vary smoothly with mass 
number. Giant resonances can therefore be described with models based on the mean-field 
approximation: the two-body nucleon-nucleon correlations are replaced by the independent 
motion of nucleons in the effective one-body potential. The nucleons themselves are the 
sources of the potential. The self-consistent model of nuclear dynamics is essential for a 
correct description of ground and excited states. Self-consistent calculations ensure that the 
same correlations that are important for the ground states also determine the dynamics of 
excited states, in this particular case giant resonances. 

The details of the time-dependent relativistic mean-field model can be found in Refs. ||[7| . 
In this section we include a short outline of the basic properties. In quantum hadrodynamics 
the nucleus is described as a system of Dirac nucleons which interact through the exchange 
of virtual mesons and photons. The model is based on the one-boson exchange description 
of the nucleon-nucleon interaction. The Lagrangian density reads 

C = i)(i 1 -d-m)^ + \{da) 2 - U{a) 

- In^or + \mlu? - I + l - m y - 

- g a i)aip - g^-y ■ uip - 9 P ^1 ■ prip - e^7 • A ^ 1 - n V ■ (1) 

The Dirac spinor x/j denotes the nucleon with mass m. m a , m^, and m p are the masses of the 
cr-meson, the w-meson, and the p- meson, and g a , g^, and g p are the corresponding coupling 
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constants for the mesons to the nucleon. U(a) denotes the nonlinear a self-interaction 

U(a) = \rn\o 2 + 1 -g 2 a 3 + \g,a\ (2) 



and ri^ u , W h ', and F^ u are field tensors. 

From the Lagrangian density the set of coupled equations of motion is derived: the Dirac 
equation for the nucleons 



id t i>i 



a f-zV - g u u) - g p rp- e— ^— —A j + f3(m + g a a) 



+9u>u + g p rp + e A 



(3) 



and the Klein-Gordon equations for the meson fields 



(d 2 t - A + ml) o = -g aPs - g 2 a 2 - g 3 a 3 (4) 

(d 2 t - A + ml) uj p = gjf, (5) 

(d 2 t -A + m 2 p )p,= g p % (6) 

(pi -A)A lt = ejJT- (7) 

In the mean-field approximation only the motion of the nucleons is quantized, the meson 
degrees of freedom are described by classical fields which are defined by the nucleon densities 
and currents. The single-particle spinors ipi (i = 1,2,..., A) form the A-particle Slater 
determinant | The nucleons move independently in the classical meson fields, i.e. 

residual two-body correlations are not included, and the many-nucleon wave function is a 
Slater determinant at all times. The sources of the fields in the Klein-Gordon equations are 
calculated in the no-sea approximation M: - the scalar density 



the isoscalar baryon current 



E^i, (8) 



f = E^7^, (9) 

i=l 



the isovector baryon current 



A 

J" = E^T^, (10) 
i=i 

- the electromagnetic current for the photon-field 

i=l 

where the summation is over all occupied states in the Slater determinant | $(£)). Negative- 
energy states do not contribute to the densities in the no-sea approximation for the stationary 
solutions, but their contribution is implicitly included in the time-dependent calculation || . 
The coupled system of equations describes the time-evolution of A nucleons in the 

effective mean-field potential. The equations are equivalent to the equation of motion for 
the one-body density operator p = p(t) 

ih^p=[h D ,p], (12) 



with an initial condition for p 



Pitin) = Pin- (13) 



hp is the single-nucleon Dirac hamiltonian defined in Eq. (|3|). Starting from the self- 
consistent solution that describes the ground-state of the nuclear system, initial conditions 
are defined to simulate excitations of giant resonances in experiments with electromagnetic 
or hadron probes. For example, the one-body proton and neutron densities can be initially 
deformed and/or given some initial velocities. The resulting dynamics can be described by 
the time-evolution of the collective variables. In coordinate space for example, these will be 
the multipole moments of the density distributions. Of course, the dynamics of collective 
variables reflects the underlying single-nucleon motion in the self-consistent potential. Since 
the Dirac hamiltonian depends on the nucleon densities and currents through the solutions 
of the Klein-Gordon equations, it is obvious that the equations of motion are nonlinear. For 
a specific choice of initial conditions the nuclear system could enter into a chaotic regime 
of motion. The collective dynamics that we describe is intrinsically classical, since it is for- 
mulated as a time-dependent initial value problem, rather than a boundary value problem. 



The single-nucleon wave functions on the other hand, satisfy the Pauli exclusion principle 
at all times, i.e. on the microscopic level of single-nucleon motion the nucleus is a quan- 
tum system. The problem is therefore how to identify and quantify chaotic dynamics in an 
ensemble of nucleons that, described as quantum objects on the microscopic level, display 
classical oscillations of collective variables. 

III. TIME-SERIES ANALYSIS 

Excitations of giant resonances result in damped harmonic/anharmonic density oscil- 
lations around the equilibrium ground-state of the nucleus. Since there are two types of 
nucleons in the nucleus, protons and neutrons, a basic distinction is made between isoscalar 
and isovector oscillations. Isoscalar motion is characterized by proton and neutron densities 
oscillating in phase. The two densities have opposite phases for isovector oscillations. In 
general, shape oscillations of the density can be represented as superposition of different mul- 
tipoles: monopole, dipole, quadrupole, etc. vibrations. Depending on the type and energy 
of the experimental probe, it is sometimes possible to selectively excite different multipoles. 
In addition to shape oscillations, spin degrees of freedom can be excited in giant resonances. 
The spin of the nucleon being naturally included in the relativistic framework, it has been 
shown in Ref. || that the time-dependent relativistic mean-field model provides a consistent 
description of the spin-multipole resonances. Another type of excitations, the Gamow- Teller 
resonances, include not only spin degrees of freedom, but also rotations in isospin space, i.e. 
a neutron is transformed into a proton or vice versa. In principle these resonances can also 
be described in the framework of relativistic mean-field models, though the actual descrip- 
tion of excitations that involve so many degrees of freedom might become very complicated. 
Therefore, in the present study we consider only the most simple situation: giant monopole 
resonances (GMR),but the results of our analysis should be valid also for shape oscillations 
of higher multipole order. These are just much more difficult to solve numerically. The 
equations of motion have to be integrated in two- or three-dimensional coordinate space, 
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and the numerical accuracy of our algorithms is simply not sufficient to obtain long time- 
series that are necessary for an analysis of nonlinear dynamics. For excitations that include 
spin and/or isospin degrees of freedom, the dynamics is more involved. Nevertheless one 
expects that the present study will also provide some insight into the nonlinear phenomena 
which occur in those more complicated excitations. 

The collective dynamical variables that characterize vibrations of a nucleus are defined 
as expectation values of single-particle operators in the time-dependent Slater determinant 
|$(t)) of occupied states: multipole moments that characterize the shape of the nucleus. 
In order to excite monopole oscillations, the spherical solution for the ground-state has to 
be initially compressed or radially expanded by transforming the radial coordinate. For 
isoscalar oscillations the monopole deformations of the proton and neutron densities have 
the same sign. To excite isovector oscillations, the initial monopole deformation parameters 
of protons and neutrons must have opposite signs. For isoscalar monopole vibrations, the 
time-dependent monopole moment is defined: 

(r 2 (t)} = i($(t)|r 2 |$(t)>, (14) 

where A is the number of nucleons. The corresponding isovector monopole moment is simply 

< rl(t) >-< r 2 n (t) > . 

Fourier transforms of the collective dynamical variable determine the frequencies of eigen- 
modes. 

In Fig. [1] we display the time-series of monopole moments which represent the isoscalar 
and isovector oscillations in 208 Pb. As in our calculation of Ref. |IJ, the NL1 effective in- 
teraction has been used for the mean-field Lagrangian (effective masses for the mesons and 
coupling constants of nucleons to meson fields). This interaction reproduces the ground- 
state properties of 208 Pb, as well as experimental data on the energies of giant resonances: 
monopole, isovector dipole and isoscalar quadrupole. However, the precise values for the 
frequencies of the eigenmodes are not crucial in our study of nonlinear dynamics. 208 Pb 
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is one of the most studied nuclei, both experimentally and theoretically. The properties 
are well known, and we have selected 208 Pb because it is a heavy spherical system, with 
relatively little fragmentation of the modes, compared to lighter or deformed nuclei. The 
experimental isoscalar GMR energy in 208 Pb is well established at 13.7 ± 0.3 MeV, and the 
isovector GMR is at 26 ±3 MeV. The isoscalar mode displays regular undamped oscillations, 
while for the isovector mode we observe strongly damped anharmonic oscillations. Of course 
the time-series alone cannot determine whether the signal of the dynamical variable displays 
characteristics of chaotic motion. In Fig. |] we show the corresponding Fourier power spectra 
in logarithmic plots versus the excitation energy E = hu. There is very little spectral frag- 
mentation in the isoscalar channel, and a single mode dominates at the excitation energy of 
^11 MeV. The Fourier spectrum of the isovector mode is strongly fragmented. However, the 
main peaks are found in the energy region 25 — 30 MeV, in agreement with the experimental 
data. The frequency of the isoscalar mode provides useful information about the underlying 
dynamics, for example the compression modulus. On the other hand, very little information 
can be inferred from the Fourier spectrum of the isovector oscillations. This observation 
would be consistent with the well known fact that nonlinear systems in the chaotic regime 
do not display any useful spectral content. In order to extract more information from the 
time-series, methods of nonlinear analysis have to be used. 

For time-series that results from linear physical processes the Fourier analysis unfolds 
the characteristic frequencies which are invariants of the dynamics, i.e. they classify the 
dynamics. For nonlinear systems the corresponding analysis is somewhat more complicated. 
In order to reconstruct the dynamics of the system from the time-series of a measured 
or calculated dynamical variable, one starts by reconstructing the phase space using time 
delays. In this procedure there are two principal quantities which have to be determined: 
the time delay and the dimension of the phase space on which the attractor unfolds. In this 
section we basically follow the prescriptions of Ref. [|1(J for the reconstruction of the phase 
space. 

The time-series in Fig. |l| have the form 

8 



x(n) = x(t + nr s ) 



n = 0,l,2,... 



(15) 



where r s is the sampling time. Since in our case the time-series is calculated by numerical 
integration of a set of partial differential equations, the sampling time can be chosen arbi- 
trarily. In order to define a coordinate system in which the structure of orbits in phase space 
can be described, time-lagged variables are used 



where T is some integer that defines the time delay. The collection of time-lagged variables 
defines a vector in <i-dimensional phase space 



In general, there is no unique prescription how to choose the optimal time lag T and the 
dimension of the space d. The time delay should be chosen in such a way that x(n + jT) 
and x(n + (j + 1)T) present two independent coordinates. If the time delay is too small, 
their numerical values will be so close to make them practically indistinguishable; if it is 
too large, the two coordinates will be completely independent of each other in a statistical 
sense, i.e. no dynamics will connect their values. Of course, one also has to avoid that the 
time delay coincides with a natural period of the system. The choice of the dimension of 
the reconstructed phase space <1e (embedding dimension), depends on the dimension of the 
attractor. cIe must be sufficiently large so that the physical properties of the attractor are 
the same when computed in time-lagged coordinates, and when computed in the physical 
coordinates, which we do not know. For example, if two points in the phase space are found 
close to each other, this should result from the underlying dynamics, and not from the small 
dimension of the phase space in which the dynamics is represented. The procedure is to 
embed the time-series in a fig-dimensional phase space. The embedding dimension has 
to be equal or larger than the minimum number of dynamical variables needed to model the 
system. 



x(n + T)= x(t + (n + T)t s ) 



(16) 



y(n) = [x(n),x(n + T),x(n + 2T),...,x(n + (d-l)T)]. 



(17) 
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One possible way to choose the time delay is to consider the linear autocorrelation 
function 

N 

[x(n + T) — x] [x(n) — x] 

Cl(T) = ^ ^ , (18) 

E [x(n) - x] 

n=l 

where 

1 N 

x = — ^x(n). (19) 

n=l 

By choosing the time lag to be the first zero of the autocorrelation function, x(n + jT) and 
x(n+(j + 1)T) become, on the average over the observation, linearly independent. The linear 
autocorrelation functions for isoscalar and isovector oscillations are shown in Fig. 0. The 
normalization is Cl(T — 0) = 1. In general, when the time-series is irregular or chaotic, 
information about its past origins is lost. This means that |Cx(T)| — > as T —>■ oo, or 
the signal is only correlated with its recent past. For the isovector mode \Cl{T)\ indeed 
displays a much more rapid decrease, as compared to isoscalar oscillations. The first zeros 
of the autocorrelation function correspond to time delays of 27 fm/c and 13 fm/c, for the 
isoscalar and isovector modes, respectively. These values could be used as time delays in the 
reconstruction of the corresponding phase spaces. This method generally produces vectors in 
phase space with components that are, on the average, linearly independent. For nonlinear 
systems a more appropriate method is to use the average mutual information. This function 
can be considered as a generalization of the linear autocorrelation to nonlinear systems, and 
it tells us how much information can be learned about a measurement at one time from 
a measurement taken at another time. For a time series x(n) and time-lag T, the average 
mutual information is defined 



N 



I(T) = J2P^(n),x(n + T)) log 2 



P(x(n), x(n + T)) 
P(x(n))P(x(n + T)) 



(20) 



n=l 

The probability distribution P(x(n)) corresponds to the frequency with which any given 
value of x(n) appears. The joint distribution P(x(n), x(n + T)) corresponds to the frequency 
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with which a unit box in the x(n) versus x(n + T) plane is occupied. The information 
functions calculated from the isoscalar and isovector time-series are displayed in Fig. [|. We 
notice that on the average there is much more mutual information in the isoscalar signal. 
The prescription is now to choose as time-lag the value for which I(T) displays the first 



minimum [II]. The curious result is that from the average mutual information we obtain 
exactly the same time delays as from the linear autocorrelation function: 27 fm/c for the 
isoscalar, and 13 fm/c for the isovector mode. 

In order to determine the embedding dimensions from the two time-series, we have used 



the method of false nearest neighbors ||12|| . For each vector ( |17D in dimension d, we define a 
nearest neighbor y NN (n). The Euclidean distance between the two vectors is 

R 2 d (n) = [x(n) - x NN (n)] 2 + ... + [x(n + (d - 1)T) - x NN (n + (d - l)T)f (21) 

Of course the two vectors are nearest neighbors if this distance is small, in a sense that we 
will define shortly. In dimension (d + 1) the distance between the two vectors becomes 

R 2 d+1 (n) = R 2 d (n) + [x(n + dT) - x NN (n + dT)f. (22) 

If the two points were nearest neighbors in dimension d, but now we find that R d+l (n) is 
large compared to R d (n), this must be due to the projection from some higher dimensional 
attractor down to dimension d. By going from dimension d to (d + 1), we have shown that 
the two points were " false nearest neighbors" . If the ratio 

\x(n + dT)-x NN (n + dT)\ 

is larger than some threshold, the assertion that the two vectors are nearest neighbors is 
false. The fact that they are found to lie close to each other in dimension d is not a 
property of the dynamics of the system, but the result of projecting the dynamics onto a 
phase space of too low dimension. The embedding dimension is now determined with a 
simple procedure. First we decide what value for the Euclidean distance should be taken 
as small. This generally depends on the data set that we are analyzing. For all vectors in 
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the phase space of dimension d we count the number of nearest neighbors. Then we have 
to determine the percentage of nearest neighbors that turn out to be false when going to 
dimension (d + 1). The minimal necessary embedding dimension dE is selected to be the 
one for which the percentage of false nearest neighbors goes to zero. For the isoscalar and 
isovector time-series, the false nearest neighbors are displayed in Fig. [5] as functions of the 
phase space dimension. The percentage of false nearest neighbors goes to zero for dE = 3 
(isoscalar mode), and for dE = 4 (isovector mode). These values are taken as embedding 
dimensions for the reconstruction of the corresponding phase spaces. 

The reconstructed phase space can be represented by the recurrence plot. By embedding 
the time-series we create a sequence of vectors 

y(n) = [x(n), x(n + (d E - 1)T)] 

in the phase space of dimension dE- In our example the time delay is 27 fm/c for the 
isoscalar, and 13 fm/c for the isovector mode. The corresponding embedding dimensions 
are = 3 and dE = 4, respectively. We can calculate the distance between any two points 
in the phase space 

5(m, n) = \y{m) — y{n)\. 

To construct the recurrence plot we choose some distance r, and ask when \y(m) —y(n)\ < r. 
m is placed on the horizontal axis, n on the vertical axis, and a dot is placed at the coordinate 
(m,n) if \y{m) — y(n)\ < r. For a periodic signal the recurrence plot displays a series of 
stripes at 45 degrees. If a time-series is chaotic, the recurrence plot has a more complicated 
structure. It is non-uniform and boxes of dense points appear along the diagonal. Of course, 
in all recurrence plots there is a stripe on the diagonal m = n. The recurrence plots for 
the phase spaces of the isoscalar and isovector time-series are shown in Fig. |6|. We notice a 
pronounced difference between the two modes. For the isoscalar mode the recurrence plot 
displays a pattern representative for regular oscillations, with stripes separated by a distance 
that corresponds to the period of oscillations. On the other hand, the recurrence plot for 
the isovector mode indicates non-stationarity. 
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The number of dots in a recurrence plot tells how many times the phase space trajectory 
came within distance r of a previous value. A measure of the density of dots is provided by 
the correlation integral. If the dynamics of a system is deterministic, the ensemble of phase 
space trajectories converges towards an invariant subset of the phase space - the attractor. 
For chaotic dynamics the attractor has fractional dimension, whereas the dimension is in- 
teger for regular dynamics. The correlation dimension of the attractor can be numerically 
evaluated from the correlation integral. If there are N points y(n) in the reconstructed phase 
space of dimension d, we can compute all distances \y(m) — y(n)\. The correlation integral 



is defined |T3| 



2 N 

° 2{r) = N(N-l) ^ @(r " W(m) ~ (24) 

for a distance r in phase space. Q(x) = if x < and Q(x) = 1 for x > 0. In a certain 
range of r, the scaling region, Cz{f) behaves like 

C 2 (r) = r d . (25) 

The correlation dimension D 2 is determined by the slope of the log C 2 (r) versus log r. It 
is defined as the slope of the plot in the r — > limit. The dimension of the attractor can 
be determined by plotting log ^(r) versus log r for a set of increasing dimensions of the 
phase space. As the embedding dimension increases, the correlation dimension D 2 should 
saturate at a value equal to the attractor's correlation dimension. The logarithm of the 
correlation integral is plotted in Fig. |7| for the isoscalar and isovector modes, for a set of 
increasing dimensions (direction of the arrow). The corresponding correlation dimensions 
are displayed in Fig. |8|. We notice that for the isoscalar mode, for d > 3, the correlation 
dimension saturates at D 2 = 2. The integer value for the dimension of the attractor indicates 
regular dynamics. For the isovector mode the correlation dimension does not saturate, but 
slowly increases to some fractional value between 2 and 3. The fractional dimension of the 
attractor would imply chaotic or stochastic dynamics. 
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IV. INFORMATION-THEORETIC FUNCTIONALS 



In this section we continue with the analysis of nonlinear dynamics in time-series of 
giant monopole oscillations. The identification and quantification of the underlying regular 
or chaotic dynamics will be based on the evaluation of information-theoretic functionals. We 
start with the information entropy of a physical system. For classical systems the entropy is 
defined on the phase space density, and represent the missing information about which fine- 
grained cell of the phase space the system occupies. For a quantum system the information 
entropy can be defined on the density operator. The von Neumann entropy of the one-body 
density operator D (measured in bits) 

S = -tr(D log 2 D) = - £ Xj log 2 Xj, (26) 

j 

where Xj are eigenvalues of the density operator (occupation probabilities), and the entropy 
can be interpreted as the missing information about which eigenvector the system occu- 
pies. Entropy is conserved under Hamiltonian dynamical evolution, both classically and 
for a quantum system. Classically this is a consequence of the conservation of phase space 
volume, in quantum mechanics it follows from the unitarity of Hamiltonian evolution. How- 
ever, if the system is coupled to a perturbing environment, the interaction generally changes 
the system's entropy. In our model only the motion of nucleons, protons and neutrons, 
is quantized. The equation of motion (|i~2"D describes the time evolution of the one-body 
density in the time-dependent meson fields. The mean-fields play the role of the environ- 
ment, the self-consistent interaction of the nucleon with the meson fields determines the 
nonlinear dynamics. This implies that, if from the nucleon densities we define some time- 
dependent entropy functionals, their time evolution might contain useful information about 
the dynamics of the system. For example, we can define the information entropy functional 

S{t) = - J p(r, t) log 2 p(f, t) d\ (27) 

where p is the vector density (0-th component in Eq. Since we consider both isoscalar 
and isovector motion, the density in ([27]) can be the neutron or the proton density, or the 
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total nucleon density. 

In Fig. |9] we display the time-dependent total (a), neutron (b), and proton (c) entropies 
( p?D for the isoscalar and isovector monopole oscillations in 208 Pb. In addition, in all three 
cases we include the entropy that results from the time evolution of the system that has not 
been excited in any way. This ground state entropies provide a measure of the numerical 
accuracy of the integration algorithm. We notice that for the isoscalar mode the entropies 
display regular oscillations which reflect the exchange of energy between the nucleons and the 
meson fields. The oscillations are identical to those of the dynamical variable, the isoscalar 
monopole moment in Fig. |lj For the isovector mode the entropies, in addition to somewhat 
more complicated oscillations, slowly decrease to the values that are characteristic for the 
ground-state of the nucleus. This decrease in entropy is caused by the strong mean-field 
damping of the isovector mode, i.e from the collisions of the nucleons with the moving wall 
of the nuclear potential generated by the self-consistent meson fields. In the isovector mode 
the protons and neutrons effectively move in two self-consistent potentials that oscillate out 
of phase, and which in this way inhibit the resonance. To extract the information content of 
the time-dependent entropies (27), we have calculated the Fourier transforms. For the total, 



neutron and proton entropies, the Fourier power spectra are displayed in Figs. ITU] and [TT 
for isoscalar and isovector oscillations, respectively. For the isoscalar mode the information 
content of the entropy is exactly the same as that of the dynamical variable, the monopole 
moment: a single mode dominates, at a frequency which corresponds to the excitation 
energy of the giant monopole resonance. This is not surprising, if one considers that the 
monopole moment is defined with an integral identical to the one that defines the entropy 
in fl2"7|), except that — log2 p is replaced by r 2 . The situation is different for the isovector 
mode (Fig. [11]). The entropy contains more information than the dynamical variable. In 
addition to the frequencies in the region of isovector monopole resonances (25 — 30 MeV), 
there are strong peaks at the frequency of the isoscalar resonance. They are related to the 
compressibility modulus of the nuclear matter. The entropy of the total density therefore 
contains information about both modes, but now we notice that the Fourier spectra for 
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the neutron and proton entropies are different. For neutrons the peaks in the region of 
isovector excitations are strongly suppressed, and there is fragmentation at the frequency of 
the isoscalar mode. 

The radius of a heavy spherical nucleus like 208 Pb is ~ 5 — 6 fm. The giant multipole 
resonances represent collective oscillations of the proton and neutron densities, and therefore 
provide excellent physical examples for the analysis of systems that have spatial as well as 
temporal structure. For a nonlinear system in chaotic regime, we might consider the influence 
of spatial motion on temporal chaos. We ask what are the spatial correlations in a finite 
system that displays chaotic oscillations of a collective dynamical variable. Consider, for 
example, the conditional entropy defined from a two-body total density 

s 2 (t) = - I P \f,f',t) log (^r'-v 



d 6 r d 6 r', (28) 
where the two-body density matrix is defined from the Slater determinant of occupied states 
P 2 = E (i\Km(k\p(r 'MmMatWiMt))- (29) 

kijl 

In coordinate representation the expression becomes 

($(f)| : P (^p(r') : |$(f)) = p{f)p{r') - 

Z N 

^^+(r)^-(r)^ + (f ')ipi(r') - E ^ti^ji^t^ 'M(0- 

(30) 

ipi{f) denotes the single nucleon Dirac spinor, and Z (N) is the number of protons (neutrons). 
The conditional entropy ( p8|) should provide a measure of two-body spatial correlations. For 
some kind of collective motion, regular or chaotic, this function contains the following infor- 
mation: how much are the oscillations of nucleon density at some point in space determined 
by the oscillations at some other point in the system, i.e. what are the correlations between 
oscillations of nucleon density at various points in the finite system. 

The time- dependent entropies ( PB[ ) that corresponds to isoscalar and isovector oscillations 
are shown in Fig. 0. They are compared with the value that results from the time-evolution 
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of the system that has not been excited (time-dependent entropy of the ground state). For 
the isoscalar mode, regular modulated oscillations are observed. Comparing also with the 
reference ground-state entropy, we notice how the numerical accuracy affects the results for 
long times of integration (T > 2000 fm/c). The entropy that corresponds to the isovec- 
tor mode is much lower and more irregular at the beginning, but it eventually approaches 
values comparable to those of the isoscalar mode. Similar to the entropy defined on the one- 
body density operator, this behavior reflects the strong mean-field damping of the isovector 
oscillations. The information contents of the "two-body" entropies are shown in the cor- 



responding Fourier power spectra in Fig. 13. For the isoscalar mode we again find that 



the entropy contains the same information as the dynamical variable, a single mode at the 
frequency of the giant resonance. This means that there is a high degree of two-body corre- 
lations for the isoscalar mode, the nucleon density oscillates with the same frequency at all 
points in the nuclear system. For the isovector mode we do not find any useful information 
in the Fourier spectrum. There is a highly fragmented structure in the region of the isoscalar 
giant resonance, but in addition we find strong peaks in the very low frequency region. This 
result indicates that there is very little spatial correlation for the isovector oscillations of the 
nucleon density, or that the nonlinear nuclear system oscillates in a regime for which the 
Fourier spectrum of the conditional entropy ( P8f) does not contain useful information. 

In the previous section, we have used the average mutual information function to deter- 
mine the time delay for the reconstruction of the phase space. This function quantifies the 
information that is contained in the signal, at some moment in time, about the value of the 
dynamical variable at other times. Since we describe isoscalar and isovector oscillations, i.e. 
we distinguish between proton and neutron components of the system, we might ask how 
much information is contained in the dynamical variable of the neutron distribution, about 
the proton subsystem, and vice versa. The two dynamical variables in this example are the 
mean square radii of the two distributions. We define the information function 

= lo £2 Pi{x^ u) ). (31) 

i 
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The signal x is quantized in units of e. The probability distribution Pi(x) corresponds to 
the frequency with which any given value of x appears in the box i of dimension e. The sum 
is over occupied boxes of dimension e, in the one-dimensional embedding of the time-series. 
For two time-series, the corresponding joint information function is 

hA e ) = - Yl p i,j( x ^ Vv) log 2 Pi,j(x n , y v ) (32) 

The joint distribution Pjj corresponds to the frequency with which a box (linear 
dimension e) in the x w versus y v plane is occupied. The average amount of information about 



the variable y that the variable x contains is quantified by the mutual information 14 



M X!y (e) = 4(e) + I y (e) - I x , y {e) (33) 

Clearly, the mutual information vanishes if Pij(x,y) = Pi(x)Pj(y), i.e. if x and y are 
statistically independent. The precise value of the mutual information will of course depend 
on the size of the box e, but one should try to find a region of values for e in which M Xty (e) 
does not vary appreciably. 

In our example of giant monopole resonances, the variable x corresponds to the mean 
square radius of the proton distribution, and y to that of the neutron distribution. The 
mutual information functions (in units of bits) are displayed in Fig. [14], for the isoscalar and 
isovector oscillations. The acceptable values for e depend on the sampling of the time-series. 
For e < 0.2 fm 2 the probability distribution functions cannot be properly determined, there 
are many empty boxes, and the calculated mutual information is not useful. For larger values 
of e the calculated mutual information changes very slightly, practically with the same slope 
for isoscalar and isovector modes. Of course, the principal result is the comparison between 
the two modes: the average amount of information that (r 2 ) of the proton density contains 
about the dynamical variable of the neutron distribution is more than a factor three larger 
for the isoscalar mode. 

Another interesting possibility is to consider the mutual information as function of the 
spatial coordinate. Instead of using as dynamical variables integrated quantities like the 



mean square radii, we can follow the time evolution of the proton and neutron densities at 
various points along the radial axis (the motion is spherically symmetric). The dynamical 
variables x and y will be the values of the proton and neutron densities at each point in 
space, and we can plot the average mutual information of the densities as function of the 
radial coordinate. The results are shown in Fig. [TJ]. In addition to 208 Pb, we also display 
the mutual information for 16 and 40 Ca. These two spherical nuclei are smaller, but have 
the advantage to contain identical numbers of protons and neutrons. In all three nuclei 
the mutual information of the proton and neutron density is much higher for the isoscalar 
mode. In fact, for 16 and 40 Ca, the mutual information for the isovector mode practically 
vanishes, and no radial dependence is observed. It is somewhat higher for 208 Pb, and with 
some modulation, slowly decreases from the center of the nucleus towards the surface. The 
isoscalar mode displays a very interesting radial behavior of the mutual information. It is 
high in the nuclear volume, but there is also a pronounced minimum at the radius that 
corresponds to the surface of the nucleus. This means that there is little correlation between 
proton and neutron densities in the surface region, they oscillate almost independently. The 
mutual information increases again beyond the ground-state radius of the nucleus, but in 
this asymptotic region the densities rapidly decrease to zero. Of course, the behavior in the 
surface region is not completely unexpected. The nucleons at the surface are less bound, and 
the effective compression modulus of the surface region is different from that in the volume of 
the nucleus. For example, the nucleus 208 Pb contains 82 protons and 126 neutrons. However, 
due to the combined effects of Coulomb repulsion between protons, and the Pauli exclusion 
principle, the protons occupy practically the same volume as the neutrons. Yet the dynamics 
for the two types of nucleons seem to be very different in the surface region. The slowly 
vibrating self-consistent potentials in which the protons and neutrons move, do not average 
on the surface in the same way as in the bulk region. It is very interesting how the details of 
the underlying nonlinear dynamics emerge in the radial behavior of the mutual information 
function. 
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V. CONCLUSIONS 



In the present work we have used the time-dependent relativistic mean-field model to 
analyze the nonlinear dynamics of giant resonances in atomic nuclei. The characteristic 
properties of these collective excitations vary smoothly with the size of the nucleus, and 
therefore a self-consistent mean-field approach provides a consistent description of nucleon 
dynamics. In particular, we have analyzed the time-series of dynamical variables that char- 
acterize the giant monopole resonances: isoscalar (proton and neutron densities oscillate in 
phase), and isovector (proton density oscillates against the neutron density). The nucleons 
move in the effective self-consistent single nucleon potentials, and the equations of motion 
describe the time evolution of the one-body density. Since the time-dependent potentials are 
calculated in a self-consistent way, the model of the nuclear system is intrinsically nonlinear, 
and chaotic motion is expected for specific initial conditions. 

From the time-series of isoscalar and isovector monopole moments of 208 Pb, we have 
reconstructed the corresponding phase spaces. The time delays have been calculated from 
the average mutual information, and the embedding dimensions determined by the method 
of false nearest neighbors. The reconstructed phase spaces have been represented by re- 
currence plots. We have found that for the isoscalar mode the recurrence plot displays a 
pattern characteristic for regular oscillations, while for the isovector mode it indicates non- 
stationarity. From the reconstructed phase-spaces we have also calculated the correlation 
integrals and the corresponding correlation dimensions. As a function of the embedding 
dimension of the phase space, the correlation dimension D2 saturates at the integer value 2 
for the isoscalar mode. On the other hand, a fractional correlation dimension is found for 
the isovector oscillations. The results confirm our conclusions from Ref. [lj] that the motion 
of the collective coordinate is regular for isoscalar oscillations, and that it becomes chaotic 
when initial conditions correspond to the isovector mode. 

The nonlinear dynamics of giant resonances has also been analyzed in the framework of 
information-theoretic functionals. For the time- dependent one-body nucleon densities, we 
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have calculated the von Neumann information entropy functionals. The Fourier analysis 
has shown that the entropy of the isoscalar mode contains the same information as the 
dynamical variable. The structure is more complicated for the isovector mode, for which 
peaks are found both in the regions of isoscalar and isovector eigenfrequencies. The spatial 
correlations have been described with a time-dependent conditional entropy defined from a 
two-body nucleon density. This function enables the study of the influence of spatial motion 
on temporal chaos. From the dynamical variables that characterize the proton and neutron 
distributions, i.e. the mean square radii, we have calculated the average mutual information 
for the isoscalar and isovector modes. The average information that is contained in the 
collective dynamical variable of the proton density, about the neutron density, is more than 
a factor three larger for the isoscalar mode. We have also analyzed the mutual information 
between proton and neutron densities as a function of the spatial coordinate. It has been 
shown that, not only is the average mutual information much higher for the isoscalar mode, 
but it also displays an interesting radial dependence which reflects the differences in the 
dynamics of the monopole motion in the volume and on the surface of the nucleus. 

The results of the present analysis, as well as those of Ref . |TJ , have shown that giant res- 
onances in nuclei provide excellent examples for the study of regular and chaotic dynamics 
in quantum systems. In addition, the finite spatial extension of nuclei enables the analysis 
of spatio-temporal behavior in nonlinear dynamical systems. And yet we have only exam- 
ined the most simple modes of collective motion: monopole oscillations. More complicated 
excitations, especially those involving spin and isospin degrees of freedom, would certainly 
disclose more interesting properties of the underlying nonlinear dynamics. 
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Figure Captions 

FIG. 1. Time-dependent isoscalar < r 2 > (a), and isovector < r 2 > — < r 2 > (b) monopole 
moments for 208 Pb. 

FIG. 2. Fourier power spectra for the isoscalar (a), and isovector (b) monopole oscillations in 

208p] 3 

FIG. 3. Linear autocorrelation functions for isoscalar (a), and isovector (b) oscillations. 

FIG. 4. Average mutual information as a function of time delay for monopole oscillations. 

FIG. 5. Percentage of false nearest neighbors as a function of phase space dimension, for isoscalar 
(a), and isovector (b) modes. 

FIG. 6. Recurrence plots for the time-series of isoscalar (a), and isovector (b) monopole oscilla- 
tions. 

FIG. 7. Sequence of correlation integrals for isoscalar (a), and isovector (b) monopole oscillations, 
log C*2(r) vs log r plots are displayed for a sequence of embedding dimensions d = 1, 8. 

FIG. 8. Correlation dimension D2 as function of embedding dimension, for isoscalar and isovector 
oscillations. 

FIG. 9. Time-dependent entropy functionals ( p7|) for isoscalar and isovector monopole motion. 
We display the total entropy of the nucleon system (a), the neutron (b), and the proton entropy 
(c). Solid curves correspond to isoscalar oscillations, dashed curves to isovector oscillations. The 
thick solid curves are the reference ground-state entropies. 

FIG. 10. Fourier power spectra of the total (a), neutron (b) and proton (c) entropies, for isoscalar 
monopole oscillations. 
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FIG. 11. Fourier power spectra of the total (a), neutron (b) and proton (c) entropies, for isovector 
monopole oscillations. 

FIG. 12. Time-dependent conditional entropies (|28| ) for isoscalar and isovector monopole motion. 
The thick solid line is the reference ground-state entropy. 

FIG. 13. Fourier power spectra of the conditional entropies (p8|) for isovector (a), and isoscalar 
(b) monopole oscillations. 

FIG. 14. Mutual information ( |33|) between the time-dependent mean square radii of the proton 
and neutron density distributions. The two curves that correspond to isoscalar and isovector 
oscillations, are plotted as functions of the size of the box in the linear embedding of the time-series. 

FIG. 15. Radial dependence of the mutual information between proton and neutron density 
distributions. Results for 16 O, 40 Ca and 208 Pb are displayed. Solid curves correspond to isoscalar 
oscillations, dashed curves to isovector oscillations. 
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